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ABSTRACT 

Aims. In this work, a sample of v sin i of B9 to F2-type main sequence single stars has been built from highly homogeneous v sin i parameters 
determined for a large sample cleansed from objects presenting the Am and Ap phenomenon as well as from all known binaries. The aim is to 
study the distributions of rotational velocities in the mass range of A-type stars for normal single objects. 

Methods. Robust statistical methods are used to rectify the vsin; distributions from the projection effect and the error distribution. The 
equatorial velocity distributions are obtained for an amount of about 1100 stars divided in six groups defined by the spectral type, under the 
assumption of randomly orientated rotational axes. 

Results. We show that late B and early A-type main-sequence stars have genuine bimodal distributions of true equatorial rotational velocities 
due probably to phenomena of angular momentum loss and redistribution the star underwent before reaching the main sequence. A striking 

I lack of slow rotators is noticed among intermediate and late A-type stars. The bimodal-like shape of their true equatorial rotational velocity 

' distributions could be due to evolutionary effects. 
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1. Introduction 

The observed rotation of stars undergoing the early main se- 
quence (MS) evolutionary stages, known as the "dwarf" stage, 
is the result of a long process of global angular momentum 
loss and angular momentum redistribution inside the stars. This 
process begins with the stellar formation phases, which encom- 
pass the fragmentation of rotating clouds (Bodenheimer 1981; 
Matsumoto et al. 1997) and the angular momentum losses re- 
quired to solve the angular momentum problem through the ini- 
tial magnetic braking (Mouschovias & Morton 1985a,b) and 
bipolar overflows (Pudritz 1985). They give the stars the char- 
acteristics that can be considered as their "initial conditions" 
at the birth-line. Further pre-main sequence (PMS) processes 
of angular momentum exchange and/or magnetic locking with 
the accretion disc (Konigl 1991; Edwards et al. 1993), angular 
momentum losses by stellar winds (Schatzman 1962; Shu et al. 
2000) and internal angular momentum redistribution phenom- 
ena, shape the internal structural patterns and rotational char- 
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acteristics of stars at the zero-age main sequence (ZAMS) that 
may be called the "initial conditions" of stars in the MS. 

Depending on the stellar internal structural pattern thus ac- 
quired, the chemical composition and the morphology of the 
magnetic fields, fossil and those created in the star (Spruit 
1999, 2002), the radiative and convective regions can be more 
or less coupled hydrodynamically (Latour et al. 1981) and/or 
magnetically (Barnes 2003) to produce the internal angular 
momentum profile (Endal & Sofia 1981; Stauffer et al. 1984; 
Pinsonneault et al. 1 990; MacGregor & Brenner 1 99 1 ; Keppens 
et al. 1995; Soderblomet al. 1993, 2001). 

Thus, according to differences in the processes and respec- 
tive efficiencies to produce loss and exchange of angular mo- 
mentum with the environment, as well as to many ways the 
stars can manage to redistribute the angular momentum in 
their interior, a variety of internal rotational laws in the ZAMS 
among objects with the same mass can be present. If only rigid 
rotation were prevalent, rather simple functions would charac- 
terize the probability distribution of a single quantity: the an- 
gular velocity rate f2/f2c (Qc is the critical rotational velocity). 
Such distributions would probably be unimodal. On the other 
hand, the conjunction of a variety of rotational laws and ro- 
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tational rates of coeval stars with the same mass in a single 
rotational distribution can have a multimodal aspect. 

Guthrie (1982) found that late B-type stars in clusters have 
bimodal rotational velocity distributions, while they are uni- 
modal for the same class of field stars. Bimodality was also 
observed for the rotation of young solar mass stars in Orion 
(Attridge & Herbst 1992; Choi & Herbst 1996; Herbst et al. 
2001; Barnes 2003). For masses intermediate to the above 
cases, Abt & Morrell (1995) found bimodal-like distributions 
among AO-FO objects in the MS, where the component owing 
to low rotators was ascribed to the chemically pecuUar Am and 
Ap stars. 

However, a high fraction of objects with the Am phe- 
nomenon were found to be binaries (Debernardi 2000). 
Regarding Ap stars, their pecuharity requires a given time af- 
ter the ZAMS to appear and thus, it can be considered some- 
how evolution-dependent (Hubrig et al. 2000; Stgpieri 2000). 
All these quoted facts can then be summarized as follows: 

1 . Although it is known that the chemical separation in Am 
and Ap stars are not produced by the slow rotation, slow 
rotation favors their appearance, so that both are correlated. 
In Am stars slow rotation can be due to their binary charac- 
ter through tidal braking, while in Ap stars could be caused 
in part by magnetic braking; 

2. Since in the whole mass interval 1.3 £ J^jJ^^ £ 3 
bimodality is clearly apparent only in cluster late B-type 
stars, it may correspond to some star formation character- 
istics. 

It is tempting then to review anew the velocity distributions in 
the main sequence phase ("dwarf" evolutionary state) of stars 
in the mentioned mass range, in order to detect possible sig- 
natures on differences in stellar formation characteristics. It is 
not excluded, however, that some particular signatures can be 
present in the rotational velocity distributions of those A-type 
stars that make the transition between objects with and without 
convectively unstable envelopes. These last can likely be due 
to stellar structure properties rather than to formation circum- 
stances. 

Thus, the aim of this paper is to see: 1) whether single late 
B- and A-type field stars in the early main sequence evolu- 
tionary phases have unimodal or multimodal rotational distri- 
butions; 2) whether the rotational distributions of intermediate 
and late A-type stars bear signatures related with the complex- 
ity of their stellar envelope/atmospheric structure. 

The data used in this paper and the selection of samples are 
described in Sect. 2. The v sin i distributions are presented in 
Sect. 3. Section 4 gives details on the statistical processing of 
stellar samples to get the equatorial velocity distributions from 
the observed v sin i values and describes the resulting distribu- 
tions. Finally, the results are discussed in Sect. 5 and suimna- 
rized in Sect. 6. 

2. Rotational velocities data 

2.1. Main v sin / data sources 

When studying distributions of stellar parameters, not only 
large samples, but also as homogeneous as possible are suited 



to avoid or minimize possible biases. The hitherto largest Ust 
of rotational velocities is provided by Glf bocki & Stawikowski 
(2000), who gathered v sin i parameters from the literature of 
over 1 1 000 stars of all spectral types without magnitude Um- 
itations. Nonetheless, this catalog is built up from about 200 
different bibhographic sources and it can be therefore highly 
heterogeneous. We preferred then to base our study on two ho- 
mogeneous data sets: on one hand, the v sin / obtained using the 
Fourier transform method (FT) by Royer et al. (2002a,b, here- 
after Papers I and n, or I U II) and on the other hand, the v sin i 
determined by Abt & Morrell (1995, hereafter called AM) 
and Abt et al. (2002, hereafter called ALG, or AMALG when 
both papers are referred to) from the FWHM of He i and Mg n 
lines calibrated into the vsin; scale of Slettebak et al. (1975). 
Nevertheless, these two data sets are different. The AMALG 
sample is complete down to the apparent magnitude V = 6.5 
mag and contains A- and B-type stars from the northern hemi- 
sphere of all classes. The I U II sample concerns objects from 
both hemispheres, some down to magnitude V = 8 mag, and 
observation of those with missing v sin ; parameters were priv- 
ileged. The selection criteria used to estabhsh the I U II stellar 
sample are suimned up in Gerbaldi & Mayor (1989). 
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Fig. 1. Distribution of spectral types of the bright stars (V < 8) with 
known v sin ( in the spectral type range B9-F2. The thin-dashed his- 
togram represents the stars from Glgbocki & Stawikowski (2000) and 
the thick-dashed histogram is the sample from Royer et al. (2002b). 
The gray histogram is for B9-type stars in Abt et al. (2002). 

The merging of v sin i data from the I U II sample and AM 
has already been dealt with in Paper II. The distribution of 
the number of stars with measured v sin ; in the resulting sam- 
ple is displayed in Fig. 1 (thick-dashed histogram). As a com- 
parison, the same distribution for Glfbocki & Stawikowski's 
in the spectral type range B9-F2 is over-plotted (thin-dashed 
histogram). In this figure we see that the exhaustive compila- 
tion by Glfbocki & Stawikowski contains a significantly higher 
amount of data for stars in the B9 to AO-type spectral interval 
and for F2-type stars. The lack noted of late B-type stars in 
I U II can, however, be nicely completed with the ALG data to 
obtain a merged sample that becomes comparable in size with 
the exhaustive collection in Glgbocki & Stawikowski (2000) in 
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Table 1. Highlight on the foiir common stars between Slettebak et al. (1975, SCBWP), Abt et al. (2002) and Papers I and n (standard deviation 

is indicated; dash "-" stands for only one measurement) and comparison with data from the literature, v sin ; from the literature are classified in 
three subgroups according to the way they are derived: by-product of a spectrum synthesis or frequency analysis of the lines profiles. Flags from 
HIPPARCOS catalog are indicated: variability flag H52 (C: constant, D: duplicity-induced variability, M: possibly micro-variable, P: periodic 
variable, -: no certain classification) and double annex flag H59 (C: component solution, -: no entry in the Double and Multiple Systems 
Annex). 



Name 


HD 


Sp. type 


SCBWP 


lun 


vsin/ (kms"') 

literature 

spec, synth. freq. analysis 


HIPPARCOS 

- H52 H59 


jSTau 


35497 


B5IV 


60 


64 ±5 






THer 


147394 


B9I1I 


30 


46 ±2 


32(1) 32(2) 


P 


yLyr 


176437 


B9.5I1I 


60 


72 ±2 




M 


0? Aqr 


222661 


B9V 


120 


150- 




C 



(i> Adehnan (1988) Smith & Dworetsky (1993) 



the whole spectral range from B9 to F2 as shown by the gray 
shading in Fig. 1. So, the rotational velocity distributions stud- 
ied in the present paper are based entirely on the I U II plus 
AMALG merged and homogenized sample of v sin i parame- 
ters. 

2.2. Merging with Abt et al. (2002) and tlie v sin / scale 

AMALG determined v sin i for nearly 2800 B- and A-type stars 
in the northern hemisphere using the full width at half maxi- 
mum (FWHM) of the He: 4471 and Mgii 4481 A Hues cali- 
brated in the v sin i scale of Slettebak et al. (1975). The merging 
with the A-type part of the sample (AM) was carried out in the 
previous paper. The intersection of the B-type part (ALG) with 
I U II contains 64 common stars, essentially of spectral types 
B8, B9 and B9.5. The comparison of both v sin i scales is dis- 
played in Fig. 2. From this diagram we can see that there is a 
systematic deviation, in the sense that the v sin i k 100 km s"' 
from I U II are on average « 16% larger than those in ALG. 
This confirms the global tendency of v sin ; being lower when 
determined with the Slettebak et al. calibration (cf. Royer et al. 
2002b). Both scales coincide, however, for v sin? ^ SOkms"^ 
The linear regression line between both scales was determined 
with GaussFit (Jefferys et al. 1998a,b), a robust least squares 
minimization program to obtain empirical functions: 

vsinj'iuii = 1. 2 1 ±0.04 V sin; ALG - 6.1±i.3. (1) 

In papers I U II there are also four standard stars common 
with Slettebak et al. which are listed in Table 1 and whose v sin i 
can be compared. Other details on these objects are given in 
Appendix A. 

It is worth emphasizing that Eq. 1 was used only for stars 
of late B spectral types, i.e. B9 and B9.5. The resulting I U II 
plus AMALG sample contains thus 2167 vsin; parameters of 
stars from spectral type B9 to F2. 

Let us comment briefly the reason of aligning the v sin / 
parameters to our I U II rotational velocity scale. According 
to the gravity darkening (GD) efl'ect, the fastest rotating re- 
gions (equatorial zone) contribute less to the global line flux 
than higher latitudes, in particular when it concerns He i 447 1 
and Mg n 4481 A lines. It has been shown (cf. Stoeckley 1968; 




V sm ijLG 

(km s"') 

Fig. 2. Comparison of v sin i data for the 64 common stars from Papers 
1 and II (1 U II) and Abt et al. (2002, ALG). The solid line stands for 
the first bisector. 



Fremat et al. 2005) that the interpretation of spectral lines of 
fast rotating stars with classical methods lead systematically to 
underestimated values of v sin j . Nevertheless, Slettebak et ill. 
(1975), who take GD into account, report that their vsin; pa- 
rameters are systematically lower than those issued from clas- 
sical methods by roughly 15 % among early stars up to spec- 
tral type A2, which is contrary to what is predicted by modern 
calculations of non-LTE rotationally broadened lines affected 
by GD (Fremat et al. 2005). The difference noted by Slettebak 
et al. (1975) drops to only some percents for stars later than 
A2. Howarth (2004) has also re-examined recently Slettebak's 
et al. V sin i scale and confirmed the above mentioned underes- 
timation. 

Since the v sin ; parameters in papers I U II are systemati- 
cally higher than those by AMALG, which mirror Slettebak's 
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et al. (1975) scale, we adopted a conservative position and de- 
cided to use our FT scale of v sin / for the merged 1 U 11 plus 
AMALG data source. Except for the assumptions that under- 
lie the FT method, the I U II v sin i scale is independent from 
any other calibration. This scale has been moreover confirmed 
recently by a new more robust FT method (Reiners & Royer 
2004) and proved to be consistent also with new v sin i deter- 
minations (e.g. Erspamer & North 2003; Fekel 2003). 

2.3. Sample definition 
2.3.1. Luminosity classes 

We attempt at studying the characteristics of the v sin i distribu- 
tions in the MS evolutionary phase. Since we cannot determine 
ages for all studied stars, we proceed to a statistical selection of 
those objects that could be meant to be in the MS stage. To this 
end, we make a selection of "dwarf" stars and retain only those 
with luminosity classes ranging from V to IV. This luminos- 
ity class criterion is large enough to take into account possible 
rotationally induced effect in fast rotators, which make stars 
look more luminous or evolved than they really are (Collins & 
Sonneborn 1977; CoUins et al. 1991; Fremat et al. 2005). The 
luminosity classes used for this selection are the ones compiled 
in the HIPPARCOS catalog (ESA 1997). 



between spectral class and B-V color index is shown in Fig. 3 . 
The mean error on B-V, &b-v - 0.0067 mag, is small in com- 
parison with the color overlap between the different spectral 
classes. A running standard deviation (on 200 points) on this 
relation enables us to estimate the error committed on the spec- 
tral classification: about one sub-class foiB-V £ 0.15 mag, 
and about two sub-classes for B - V ^ 0.15 mag. The color 
index B - V - 0.15 mag approximately corresponds to the 
spectral class A5 (Cowley et al. 1970). Color indices can be 
rather strongly affected in late type stars by fast rotation effects 
if rotational rates are Q/Qc ~ 0.9 (Collins & Sonneborn 1977; 
Collins et al. 1991; Fremat et al. 2005). Since we cannot cor- 
rect in advance individual spectral types for possible induced 
rotational changes, we shall conclude a posteriori that this phe- 
nomenon may in principle be aff'ecting only a low fraction of 
stars (see Fig. 10: 7.6 % in B9, 5.4 % in AO-Al, 4.2 % in A2- 
A3, 4.2% in A4-A6, 3.6% in A7-A9 and 2.5% in F0-F2). 
The full spectral range B9-F2 was thus divided into six groups, 
which are listed in Table 2. 

Table 2. Subsamples of Main-Sequence stars (classes V, IV-V and IV): 
the mean radii R/Rq and the mean masses ^/^q are derived from 
Habets & Heintze (1981), as well as the relative variation of the radius 
AR and of the mass along the interval of spectral type. For each 
subsample, the number of stars is given. 



2.3.2. Spectral types 
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Fig. 3. Comparison between the spectral class and the B - V color in- 
dex (both from HIPPARCOS catalog) for the stars in the v sin ; sample, 
with luminosity classes V, IV-V and IV. The colored contour indicates 
the number of stars per bin of 0.01 mag, per spectral class: black — 1 
to 5 stars; dark gray — 5 to 10 stars; light gray — 10 to 15 stars; white 
— more than 15 stars. The average error bar in B - V is indicated: 
as-v = 0.0067 mag. 

The selected 1500 stars were finally divided into spectral 
type subsamples, based on the spectral types compiled in the 
HIPPARCOS catalog. In establishing the subgroups, we looked 
for a compromise between the sample size and the radii depen- 
dence with mass. The correspondence between spectral type 
and radius used is from Habets & Heintze (1981). The relation 



Spec, types R/Rq 


AR(%) 




(%) # stars 


B9 


1.95 




2.34 


198 


AO-Al 


1.79 


±0.8 


2.04 


±1.5 397 


A2-A3 


1.70 


±1.2 


1.90 


±2.1 377 


A4-A6 


1.60 


±2.5 


1.73 


±3.8 182 


A7-A9 


1.48 


±2.7 


1.52 


±4.3 190 


FO-1-2 


1.37 


±2.2 


1.37 


±2.2 197 



2.3.3. Chemically peculiar and binary stars 

In order to clean the sample from stars whose rotation could 
have been modified by tidal or magnetic braking, chemically 
peculiar stars (CP) and "close" binary stars (CB) have been 
separated from "normal" stars. This choice is moreover jus- 
tified because of the a priori selection biases in the I U II sam- 
ple. As regards the southern sample from Paper 1, selection 
criteria are described by Gerbaldi & Mayor (1989). They fo- 
cused on "normal" A-type stars, excluding any known peculiar 
or spectroscopic binary stars. The same criteria have been used 
for the northern sample in Paper II. 

CP stars (chemically peculiar stars): The catalog of Ap and 
Am stars from Renson et al. (1991) is used to identify pecu- 
liar stars, as well as the spectral classification made by AM and 
ALG. The number of CP stars in each subsample is summa- 
rized in Table 3. 



CB stars ("close" binary stars): The tidal efi'ect in a multiple 
system tends to synchronize the axial rotation period on the 
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Fig. 4. Distributions of V magnitudes and completeness of the CP 

and CB stars. The numbers of stars are given in logarithmic scale. 
Top panel) the distributions (in the range 3-8 mag) are displayed for 
the A-type stars in the HIPPARCOS catalog, and their intersection 
with the lists given by Renson et al. (1991) and Pedoussaut et al. 
(1985). Middle panel) the distributions are shown for the vsin/ sam- 
ple (hatched histogram) and the selected CP stars (light gray). Dark 
gray histogram is the v sin i sample scaled to the number of CP stars in 
the V-magnitude range 4-6 mag. Bottom panel) the distributions are 
shown for the v sin i sample (hatched histogram) and the selected CB 
stars from the list by Pedoussaut et al. (1985) only (light gray). Dark 
gray histogram is the v sin / sample scaled to the number of binary 
stars in the V-magnitude range 4-5 mag. 



orbital period. This category of stars was selected on criteria 
based on HIPPARCOS data and spectroscopic data. Except for 
a few of them, all selected stars are in the HIPPARCOS catalog 
(ESA 1997). The binaries detected by the satellite with AHp < 
4 mag are flagged as CB stars. The compilation by Pedoussaut 
et al. (1985) of spectroscopic binaries is used to complete the 
identification, together with the "Eighth catalog of the orbital 
elements of Spectroscopic Binaries" (Batten et al. 1989). 



In what follows, all stars that do not obey to CB or CP cri- 
teria are simply called "normal". Since these "normal" stars 
constitute the bulk of our stellar sample, it is worth having a 
statistical inference of the fraction of stars where the CP and 
CB characters are stiU unknown. 

Completeness of the knowledge of CP stars: a large part of 
the data sample was observed by Abt (AMALG), who derived 
at the same time the v sin / values, determined the spectral clas- 
sification and identified the chemical pecuUarities. In addition, 
the catalog of Renson et al. (1991) was also used to account 
for peculiar stars. As a reference to test the completeness of 
this catalog we plot in the top panel of Fig. 4 all A-type stars 
in the HIPPARCOS catalog, the CP stars in the Renson et al. 
(1991) catalog and Pedoussaut's et al. (1985) compilation of 
spectroscopic binaries. The number of stars is displayed in log- 
arithmic scale against the V-magnitude. In the middle panel of 
Fig. 4 there is the merged I U II plus AMALG data set (hatched 
histogram), where the short-dashed line is simply a fit of their 
distribution up to V = 6.5 mag. The light gray histogram gives 
the CP stars identified in this paper over which we have super- 
imposed the long-dashed Unear fit valid up to V = 6.5 mag. We 
see that the fraction of CP stars to all stars in the merged sam- 
ple is fairly constant for all magnitudes brighter than V - 6.5 
mag and it represents roughly 15 % of objects. The same frac- 
tion, and, in the same magnitude interval, is also seen in the top 
panel. Making the hypothesis that the proportion of CP stars in 
our merged ensemble remains at 15 % also for V ^ 6.5 mag, the 
scaled distribution (dark gray) shows the deficit of undetected 
CP, which amounts only to some 20 stars. 

Completeness of the knowledge of CB stars: the identifi- 
cation of the binary nature of stars is a difiicult problem and 
can lead to partially satisfying solutions. Most of the time it 
requires surveys with spectroscopic monitoring. The use of 
HIPPARCOS results to account for close binary stars sim- 
ply does not allow detection of the closest ones (p < (X'l- 
(y.'15, ESA 1997) or those with large magnitude difference 
(AHp > 3. 5^ mag). Considering the mass-luminosity relation 
from Gerbaldi et al. (1999), the magnitude difference leads ap- 
proximately to A//p =s -2.51A(log^) =s -11.42A(log^) 
and then for AHp = 4 mag, the corresponding mass ratio is 
q ^ 0.45. 

Pedoussaut et al. (1985) provide an "exhaustive" list of 
spectroscopic binaries and the estimation of the lack of detected 
binaries is based on this list only. As it is obvious, we aim at re- 
moving binaries from the v sin ; sample to minimize rotational 
changes induced by tidal braking and the consequent synchro- 
nization. The completeness limit magnitude of the list of spec- 
troscopic binaries is estimated in the top panel in Fig. 4, which 
hardly reaches V = 5 mag. The distribution of the v sin i sample 
is scaled as in the previous paragraph, according to the number 
of spectroscopic binaries in the magnitude range 4-5 mag. The 
scaled distribution (dark gray, bottom panel) gives an estimate 
of possible undetected spectroscopic binaries, which amounts 
to some 200 stars, mainly for V ^ 5 mag. It is however worth 
noting that the frequency of binaries decreases with decreasing 
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mass ratio q (Wolff 1978) and that correspondingly the syn- 
chronization time increases (Zahn 1975). On the other hand, 
the maximum of the period distribution of spectroscopic bi- 
naries in the spectral range A0-A9 is for P « 5 d and that for 
P a! 1 d the frequency slumps to 25 % (Kogure 1981). The most 
substantial fraction of stars in the A0-A9 range with v sin i that 
can be expected to be affected by tidal braking must then cor- 
respond to q v 0.5 and P » 2 d. The synchronization time is 
a strong function of the DjR ratio, where D is the separation 
of components and R the radius of the primary: fsyn °^ (D/R)^-^ 
(Zahn 1975). For q ^ 0.5, P 2 d, average stellar mass 1.8 
and radius 1.6 Rq we obtain D/R ^ 6, so that /syn ~ 0.5 ?ms. 
where ?ms is the time spent in the main sequence. At shorter 
periods the number of spectroscopic stars decreases, while for 
longer periods the synchronization time increases strongly, i.e. 
for P a; 5 d, it is ?syn ^ 100 ?MS* This means that the possibly 
1 3 % unidentified binaries cannot blur sensitively our statistics. 
Moreover, it is easy to show that a slowing down of rotation by 
50 % of the rotational velocities in the missing binaries would 
alTccl the overall average estimate of vsin? by less than 4%, 
because more than half have P ^ 5 d. 

It is worth noting that these deficits of stars are not account- 
able only to slow rotators. Among the 218 identified "CP" stars, 
only 87 have vsin; < 70kms"', i.e. about 40%. As far as 
the spectroscopic binaries are concerned, among the 79 stars 
in common with the list by Pedoussaut et al. (1985), 43 stars 
have V sin ; < 70 km s"' , i.e. about 54 %. 

The number of objects in the different spectral-type sub- 
groups, and the different categories ("normal", CB, CP) are 
given in Table 3. All the objects are listed in Table 4. 

3. Projected rotational velocities distributions 

The histograms of v sin i of each of the six subsamples are dis- 
played in Fig. 5. The associated statistical estimators are listed 
in Table 3. Distinction has been made between the "total" sub- 
sample ans its components: "normal", "CP" and "CB" stars (re- 
spectively represented as dark hatched, light hatched and plain 
white areas in Fig. 5). The distinction in Table 3 is also made 
according to the luminosity class: "aU" classes (V, IV-V and 
IV), "V" only and "IV" only. Mean, median and dispersion es- 
timators are also given for each group. 

Several trends are noticeable in the v sin i distributions, as 
seen in Fig. 5: 

• The V sin i distributions for most massive stars (B9 and AO- 
Al) exhibit hints of a bimodality, even when only taking 
into account "normal" stars. A double peak in the v sin ; of 
AOV-type stars was pointed out by Ramella et al. (1989). 
The distributions for B9 and AO-Al groups are very sim- 
ilar. This similarity becomes quite clear if the histograms 
normalized to the number of stars in each group are com- 
pared (as they are displayed in Fig. 7). There are two 
maxima, one near vsin; ^ 160 km s"^ and the other near 
vsini ^ 50kms"^ The excess of stars with small vsin? 



Table 3. Moments of the v sin / distributions for each spectral subsam- 
ple. "total" corresponds to the whole subsample (including CP and 
CB stars); "normal", CP and CB stars are also distinguished. For both 
"total" and "normal" subsamples, the moments are given for different 
groups of luminosity class: all (IV, IV-V, V), V and IV (class IV-V 
luminosity stars only are not considered because they represent 
a very small fraction of the sample). The number of stars (#), the 
center and the dispersion of the distribution are given. Centers are es- 
timated using the median and the mean, and the dispersion is evaluated 
with the standard deviation. 



Subsample 






Center (km s ' ) 


Disp. 










median 


mean 








all 


198 


137 ±8.8 


145 ±7.0 


no . 
98 ±4.9 




total 


V 


163 


153 ±9.6 


154 ±7.7 


98 ±5.4 






1 V 




68 ±23.9 


102 ±19.1 


yj ±13.5 


B9 




all 


125 


159 ±11.0 


161 ±8.7 


98 ±6.2 




normal 


V 


1 ns 

iUo 


165 ±11.8 


167 ±9.4 


yo ±6.7 






IV 


1 1 
1 1 


103 ±32.0 


110 ±25.6 


oc .... 
03 ±18.1 




CP 




Jo 


49 ±20.4 


97 ±16.3 


98 ±11.5 




CB 




J 1 


138 ±17.6 


136 ±14.1 


86 ±10.0 






all 


397 


123 ±5.2 


129 ±4.1 


82 ±2.9 




total 


V 


363 


123 ±5.4 


131 ±4.3 


83 ±3.1 






IV 




80 ±16.2 


88 ±12.9 


OJ ±9.1 


AO-Al 




all 


271 


128 ±5.9 


131 ±4.7 


78 ±34 




normal 


V 




128 ±6.2 


131 ±5.0 


/ o ±3.5 






IV 


1 A 

14 


100 ±214 


105 ±17.1 


64 ±12.1 




CP 




oy 


123 ±14.5 


139 ±11.6 


96 ±8.2 




CB 




D 1 


107 ±13.9 


113 ±11.1 


±7.8 






all 


311 


123 ±4.5 


125 ±3.6 


69 ±2.5 




total 


V 


307 


131 ±4.9 


133 ±3.9 


68 ±2.7 






1 V 


ou 


70 ±9.0 


80 =7.2 


JO ±.1.1 


A2-A3 




all 


258 


128 ±5.2 


133 ±4.2 


67 ±3.0 




normal 


V 


ziu 


133 ±5.7 


140 ±4.5 


66 ±3.2 






IV 


A 1 

41 


71 ±11.7 


93 ±9.3 


DO ±6.6 




CP 




oU 


84 ±12.5 


105 ±10.0 


/ / ±7.1 




CB 




^0 

jy 


102 ±10.3 


109 ±8.2 


63 ±5.8 






all 


182 


133 ±5.5 


141 ±4.4 


59 ±3.1 




total 


V 


132 


158 ±6.5 


152 ±5.2 


60 ±3.7 






i V 


DO 


107 ±10.5 


107 ±8.4 


jU ±5.9 


A4-A6 




all 


137 


138 ±6.0 


144 ±4.8 


56 ±3.4 




normal 


V 

V 


07 


159 ±7.2 


157 ±5.7 


56 ±4.0 






rv 


27 


107 ±11.3 


110 ±9.0 


47 ±6.4 




CP 




19 


107 ±19.8 


106 ±15.8 


69 ±11.2 




CB 




26 


141 ±15.4 


147 ±12.3 


63 ±8.7 






all 


190 


134 ±5.5 


135 ±4,4 


60 ±3.1 




total 


V 


123 


149 ±6.5 


144 ±5.2 


57 ±3.7 






IV 


46 


110 ±11.1 


112 ±8.9 


60 ±6.3 


A7-A9 




all 


141 


133 ±6.3 


134 ±5.0 


60 ±3.6 




normal 


V 


90 


149 ±7.6 


143 ±6.1 


58 ±4.3 






rv 


32 


110 ±13.3 


114 ±10.6 


60 ±7.5 




CP 




16 


104 ±23.9 


113 ±19.1 


76 ±13.5 




CB 




33 


149 ±10.7 


149 ±8.5 


49 ±6.0 






all 


197 


107 ±5.2 


113 ±4.2 


59 ±3.0 




total 


V 


125 


117 ±6.7 


122 ±5.3 


59 ±3.8 






IV 


55 


99 ±9.7 


102 ±7.7 


57 ±5.5 


F0-F2 




all 


150 


1 10 ±5.9 


1 14 ±4.7 


57 ±3.3 




normal 


V 


91 


1 17 ±74 


123 ±5.9 


57 ±4.2 






IV 


44 


100 ±11.3 


108 ±9.0 


60 ±6.4 




CP 




18 


84 ±18.7 


92 ±14.9 


63 ±10.5 




CB 




29 


123 ±14.3 


120 ±11.4 


61 ±8.1 
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Fig. 5. Cumulative distributions of v sin i of the six subsamples of main sequence stars from B9- to F2-type. Dark hatched, light hatched and 
plain white areas, from bottom to top, respectively stand for "normal", CP and CB stars. 



Table 4. (extract) List of the 1541 B9- to F2-type stars, with their 
V sin I value, specttal type, associated subgroup and classification (CP, 
CB, blank stands for "normal"). 



HD 


V sin I 
(kms^') 


Sp. type 


subgroup 


classification 


3 


228 


AlVn 


AO-Al 




203 


170 


F21V 


F0-F2 




256 


241 


A21VA' 


A2-A3 




319 


59 


AlV 


AO-Al 


CP 


431 


97 


A7IV 


A7-A9 


CB 



cannot be considered due to Am and Ap stars, as most of 
them were taken off from the sample. 



There is a regular decrease of the maximal values 
of vsinj from B9 to F0-F2 groups, which is only 
marginally correlated with the critical equatorial veloc- 
ity Vc oc i^lR)^!^; i.e. Vc(B9)/Vc(F0-F2) ~ 1.1, while 
V sin !max(B9)/v sin J„,ax(F0-F2) ~ 1 .4. 
The proportion of low vsin/ stars 70km s"'), in the 
full subsamples, diminishes from B9 to A6-type stars (B9: 
31%; AO-Al: 31%; A2-A3: 24%; A4-A6: 10%) and 
then increases towards the later types (A7-A9: 14%; FO- 
F2: 24 %). There is a striking lack of stars with v sin ; < 
50 km s ' in the A4-A6 group. The same tendency is found, 
considering "normal" stars, but the proportion is lower. The 
effect of the proportion of low v sin i stars is also reflected 
in the variation of the dispersion (see Table 3). It decreases 
from B9 to A4 stars, and reaches a plateau around 60 km s~' 
as soon as the distribution is not significantly skewed. 
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From this point, the assumption that stellar rotation axes are 
randomly oriented is adopted. This hypothesis has been tested 
many times (Gray 1992; Gaige 1993) and is still the most vahd. 
The following sections focus on the "normal" stars, exclud- 
ing the stars classified as CP or CB in the previous section. 

4. Distributions of equatorial rotational velocities 

4.1. Rectified distributions 

The observed v sin / parameter is the projection of the equato- 
rial velocity v of the star on the line of sight, i being the in- 
clination between the stellar rotation axis and the line of sight. 
The Probability Density Function (hereafter PDF) of the v sin ; 
is thus the result of the convolution between the distribution 
of "true" equatorial velocities v, the distribution of incUnation 
angles /, and the observational error law. 

Let J? be the true projected rotational velocity of a star, the 
PDF of V sin i can be written as: 



<b(v sin i) 



I 



*?(»?) P(v smi\&)d&. 



(2) 



where P(v sin /I??) d?? is the conditional probability of v sin ; e 
+ d&]. is the PDF of "true" projected velocities, 
which itself can be written as: 



I 



(3) 



where P(§\v) dv is the conditional probability of§€ [v,v + dv], 
and T(v) is the PDF of the "true" equatorial velocities. 

The aim of this section is to recover the distribution T(v) 
from the observations. 

4.1 .1 . Smoothing of the distributions 

The observed PDF (i>(v sin /) can be estimated from the v sin / 
data. For that purpose, the kernel estimator O, with kernel K, is 
used. Details can be found in Appendix C.l. The estimator O 
is defined as: 



1 ■A 1 /x-Xi\ 

1=1 



(4) 



where n is the sample size and h is the window width, also 
called the smoothing parameter. The error on v sin i is expected 
to be multipUcative, as shown in Papers I and 11. This can 
be grasped by considering the v sin / as a scale/dilatation pa- 
rameter, like the dispersion (whereas the radial velocity is a 
shift piirameter like the mean). The error on the empirical dis- 
persion is asymptotically proportional to the dispersion (see 
Appendix B). The natural error law would thus be a log-normal 
distribution for the v sin i. In Paper I, the distribution of ob- 
served V sin ; for Sirius is perfectly fitted either by a normal or 
a log-normal distribution. This noticing triggers the processing 
of this step in logarithmic space of velocities, to get rid of the 
variance heterogeneity (Zar 1999). The data are significantly 
higher than zero to safely use ln(v sin i) instead of ln(v sin 
(Bartlett 1947). The smoothing of the PDF is then performed 



on logarithmic velocities, and the PDF Oiog(ln(v sin i)) is com- 
puted by the Kernel method, using a Gaussian kernel (standard 
normal distribution): 



m = 



1 



(5) 



The smoothing parameter h in Eq. 4 is estimated using the 
scheme described by Sheather & Jones (1991), and based on 
each subsample data. The estimator h, as solution of the Eq. 12 
in Sheather & Jones (1991), is given in Table 5 for each sub- 
sample. 



Table 5. Estimated bandwidth for the kernel method and mean inte- 
grated squared error: for each subsample, the size n (number of stars) 
is given, as well as the estimated h according to Sheather & Jones 
(1991), computed in the logarithmic velocity scale, and the variability 
band width e (Eq. C.4). 



Subsample 


n 


h 


£ 


B9 


125 


0.223 


0.0503 


AO-Al 


271 


0.147 


0.0421 


A2-A3 


258 


0.170 


0.0401 


A4-A6 


137 


0.152 


0.0583 


A7-A9 


141 


0.137 


0.0604 


F0-F2 


150 


0.149 


0.0561 



It must be emphasized, however, that using a unique value h 
for a given sample of logarithmic v sin / data will lead to create 
artificial and non significant distribution modes at low v sin i. 
For low V sin i values, the bandwidth h constant in logarithmic 
space is underestimated, due to the limitation of the precision 
on derived v sin i. In addition, for some of the subsamples, the 
number of slow rotators is small, which contributes to the pres- 
ence of non significant peaks in the distributions. These peaks 
are clearly visible in the linear velocity space (Fig. 7). 

Bowman & Azzalini (1997, Chap 2.3) propose a way of 
quantifying the variability of the density estimates and there- 
fore assessing the significance of the present modes. Using 
this method, only variance is assessed to compute the result- 
ing "variability bands", but the graphical display of this vari- 
ance structure allows the disentanglement between the genuine 
features and the "noisy" ones. The variabiUty bands are esti- 
mated by the e parameter (Table 5, see Appendix C.l) and are 
over-plotted as gray areas in Figs 6 and 7. 

The PDFs 3)iog(ln(v sin i)), for each subsample, are the re- 
sult of the kernel density estimations. The result for the largest 
subsample (AO-Al) is shown in Fig. 6. 

4.1 .2. Rectification of the error law effect 

This operation is an intermediate step in order to recover the 
PDF of true equatorial velocities. The rotational velocities and 
its associated error can be written as: 



vsinj'iavsinj « v sin j 7(1+ a), 



(6) 
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V sin i (km 
10 



100 



0.6 



cu 



0.4 - 



0.2 



0.0 




3 4 
\n{v sin i) 



Fig. 6. Distributions of the logarithmic rotational velocities ln(v sin i) 
for the AO-Al subsample. The corresponding scale of v sin (' is plotted 
as the upper x-axis. The histogram represents the data for "normal" 
stars, displayed in Fig. 5. The solid line is the smoothed distribution 
<l)iog(ln(v sin i)), obtained by using a Kernel method with a Gaussian 
kernel, on the ln(v sin i). The thick dotted line stands for the distti- 
bution yiog(lni?) of the "true" logarithmic velocities, rectified from a 
Gaussian error law (see text). 



for small values of a (i.e. terms of second order are negligible), 
the error can be written on a multiplicative form^ Therefore 
the logarithm scale can be considered: ln(vsinO ± ln(l + a). 
In this case, the coefficient a is taken to be 0.1. In Papers I 
and II, it was respectively 0.06 and 0.05 for v sin i derived from 
ECHELEC and AURELIE data, but for the v sin / measured by 
Abt, the relative error is closer to 10%. Using a - 0.1, the 
second order is obviously negligible, and a 1 %-effect on the 
observed v sin / is really not significant. 

The error of ln(v sin i) is assumed to be normally distributed 
with a mean of and a variance of (ln(l + or))^. The conditional 
probability P(ln(vsinj)|hi??) takes the shape of a Gaussian 
function: 



/•(MysinOlhii^) = 



1 



_ 1 /' In(vsinO-liu? \^ 
-e 2V InCl+or) ) 



(7) 



^|2^i: hi(l + a) 
in the equation: 

<I>iog(ln(v sin 0) = ^ 'I'logCln &) P{\n{v sin 01 hi j?) d hi &. (8) 

1'iog(lnj?) is derived using the Lucy (1974) iterative tech- 
nique. Further details are given in Appendix C.2. The PDF 
<l)iog(ln(vsin/)) and *l'iog(ln)?) of the subsample AO-Al are 
shown in Fig. 6 and illustrate the steepening and sharpening 
of the distribution deconvolved by the error law. 

The resulting distributions are transformed back in the Un- 
ear velocity space following: 



^logQni^) 



(9) 



' 7 is the symbol times/divide, multiplicative equivalent of 
plus/minus ± 



The distributions of "true" projected velocities ^(ff) are 
displayed for all the subsamples in Fig. 7. 

4.1 .3. Rectification of tine projection effect 

Under the assumption of randomly oriented rotation axes, the 
conditional probability P{§\v) in Eq 3 is: 



V Vv2 - 



(10) 



where H is the Heaviside function^. 

Combining Eqs 3 and 10, the PDF of & becomes: 



.^2 



:dv. 



(11) 



which is an Abelian integral that can be analytically solved. 
Details about this Abelian form of the convolution can be found 
in Chandrasekhar & Munch (1950). Jorissen et al. (2001) also 
describe it in the purpose of the deconvolution of the distribu- 
tion of projected masses. 

Therefore two ways are possible to compute the final dis- 
tribution of equatorial velocities T(v). In this work, only the 
Lucy-iteration technique (Lucy 1974) has been used. 

These distributions T(v) are displayed in Fig. 7 for the six 
subsamples, and described in the next Section. 

4.2. Description of tlie v distributions 

The distributions of true equatorial velocities imply the follow- 
ing facts: 

• The bimodal character of distributions corresponding to B9 
and AO-Al field dwarf stars is clearly present, even the 
samples have been cleared from all known Am, Ap and 
close binaries. For B9-type stars: the mode of slow rotators 
is ~ 35kms and for fast rotators ~ 190 kms '. For the 
AO-Al subsample: they are ~ 60km s"' and ~ 200km s"' 
respectively. The modes of fast rotators (v ^ 150 km s~^), 
for B9 and AO-Al type stars, are fitted by Maxwellian 
PDFs (Fig. 8). As shown by Deutsch (1967), the distribu- 
tions of V are well fitted by MaxwelUan functions: 



Tfast(v) 



12 In , V 
■ V e 2<. . 



j3/2 



(12) 



The parameter a is defined by the mode of the distribution 
V2fl. The proportion of slow rotators (v ^ 150kms"') is 
taken as the excess compared to the Maxwellian fast rotator 
distribution. The results are displayed in Table 6. The slow 
rotator peak corresponds to some 19 stars forB9-type stars, 
and about 66 stars for AO-Al -type stars. 
In the A2 to A9 spectral type groups, the small wiggles 
of distributions in the velocity interval ^ v ^ 70 kms"' 



^ The Heaviside step function H is defined by: 

„ , , / 1 if ;c > 
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0.008 - 



^ 0.006 - 



B9 




0.000 



AO-Al 



0.008 - 



^ 0.006 



A2-A3 




A7-A9 



200 
rotational 



300 400 500 
velocity (km s ') 




A4-A6 



F0-F2 



100 200 300 400 500 
rotational velocity (km s ') 



Fig. 7. Distributions of rotational velocities for the observed subsamples: shaded histograms are the observed v sin ; and are identical to the 
histograms for "normal" stars in Fig. 5; dashed lines stand for the smoothed distributions of true projected rotational velocities ^(ff); solid 
thick lines are the distributions of true equatorial velocities T(v) and the gray strips are their associated variability bands. Dotted lines are 
the distributions found by AM: in our AO-Al sample both peculiar and normal star distributions are overplotted, in our A2-A3 sample the 
distribution of A3-A4 stars found by Abt & Morrell (1995) is overplotted, and in our A4-A6 and A7-A9 samples the distribution of A5-F0 
stars is shown. 



Table 6. Characteristics of the bimodal distributions: the fast rotator 
mode yfla from the Maxwellian distribution parameter a, the propor- 
tion p of slow rotators as the excess compared to the Maxwellian fast 
rotator distribution, the corresponding number of slow rotators in 
the observed sample, and the mean of the slow rotator distribution /Xs. 



Spec, types 




P 


Afs 


Us 




(kms-i) 






(kms-i) 


B9 


190 


0.15 


19 


46 


AO-Al 


200 


0.24 


66 


59 



concern a negligible fraction of stars. Moreover the vari- 
ability bands associated to the distributions argue in the 



sense that the presence of these slow rotators is not signif- 
icant. These objects are probably unknown synchronized 
binaries or chemically peculiar stars that pollute the sam- 
ple of "normal" stars. There is a net lack of rotators with 
V $ 70kms An absolute minimum of this fraction is 
seen in the A4-A6 group. 
• An excess of slow rotators with velocities Q % v % 
60kms~' appears in the F0-F2 subgroup. This mode of 
slow rotators, around 20kms"', in the distribution for 
early-F stars, is significant. The B-V distribution for all the 
"normal" F0-F2 stars, ranging from 0.2 to 0.4 mag (Fig. 9), 
shows that all (but one) stars with vsin/ < 75kms"' have 
color index B-V ^03 mag. It ensues that the slow rotator 
peak {Q % V ^ 60kms"^) is only composed of stars red- 



F. Royer et al.: Rotational velocities of A-type stars. III. 



11 




V (km s ') 

Fig. 8. Slow rotators in B9 and AO-Al type stars: the PDF of equa- 
torial velocities v (solid line) are displayed. The fast rotator mode is 
fitted by a Maxwellian distribution (dashed line), and the slow rota- 
tor distribution is estimated as the excess compared to the Maxwellian 
distribution (hatched area). The difference between both distributions 
is plotted (bold solid line) together with its propagated variability band 
(gray strip). For the AO-Al type stars, the distribution of chemically 
peculiar objects found by AM is overplotted (dotted Une), normalized 
to represent the same number of stars. 



250 F 




0.20 0.25 0.30 0.35 0.40 

B-V 



Fig. 9. Projected rotational velocities as a function of the B-V index 
(From HIPPARCOS catalog) of the 150 "normal" stars in the F0-F2 
subgroup. Filled circles (•) stand for FO-type stars, open squares (□) 
for Fl and asterisks (*) for F2-type stars. The hatched strip emphasizes 
the lack of object in the area (v sin ; < 75 km s'^,B-V < 0.295 mag), 
the only object in this area is HD 36876. The dotted line is the 40-point 
running average of v sin i for the 150 objects. 



der than B - V - 0.3 mag. The drop of rotation along the 
main sequence through the early F-type stars was pointed 
out by Wilson (1966) and Kraft (1967). The decrease in ro- 
tational velocity observed in Fig. 9 is the beginning of this 
rotational break. A similar decrease occurs also for class III 
giants (Gray 1989). 



• In the high-velocity side of distributions (Fig. 7) a well 

defined peak can be noticed at v ^ 200 km s ' in all 
groups, except for the F0-F2 for which at v ^ 200 km s"' 
there is only an inflection. 

• A Roche representation of the rotating stellar surfaces is 
used to transform the A2 to F2 velocity distributions into 
Q/Qc-probabiUty distributions (see Appendix D). They are 
shown in Fig. 10. It can be seen that in each group, more 
than half of the stars rotate with Q/Qc k 0.6 and that the 
marked lacks of stars are for Q/Qc ^ 0.3. Objects with 
Q/Qc ^ 0.6 can undergo geometrical deformations induced 
by the rotation, which implies that their spectral character- 
istics may be somewhat aspect angle dependent. 

4.2.1 . Comparison with AM 

Figure 7 compares the distributions of true rotational velocities 
from AM (dotted lines) with those obtained in this work (solid 
lines). For the later A-type subgroups (A2-A3, A4-A6, A7- 
A9), there is good agreement between the distribution of nor- 
mal stars in AM with ours. It suggests that our selection of non 
CP stars is efficient. As regards our distributions for A4-A6 
and A7-A9 subgroup, compared to the corresponding A5-F0 
in AM, the tails towards lower v is much steeper in our study. 

On the other hand a strong discrepancy is observed in the 
AO-Al subgroup. The fast rotator mode found in the present 
work corresponds to the distribution found by AM for normal 
stars (Fig. 7). The excess of slow rotators emphasized in Fig. 8 
has a distribution similar to the Ap+Am stars found by AM, 
even though to discard the CP and CB stars, the same selection 
criteria as for later type stars were applied. 

Abt (2000) argues that rotational velocity alone discrimi- 
nates the normal A-type stars from the peculiar ones. However 
Mathys (2004) points out that Abt's definition of peculiar stars 
encompasses stars that do not exhibit the abundances anoma- 
lies of "classical" CP stars. 

5. Discussion 

The main objective of this paper was to obtain distributions of 
surface rotational velocities of a sample of field stars in a spec- 
tral type range where objects with the earlier spectral types 
have radiative envelopes and those with cooler spectral types 
have convective envelopes. The present research is intended 
to detect signatures of multimodal rotational velocity distribu- 
tions among "normal" dwarf stars, which are in the early MS 
evolutionary phases. 

Up to now, for stars with masses 13 ^ M jMco ^ 3.0, 
there are indications for bimodality only for cluster late B-type 
stars (Guthrie 1982) and for near solar mass stars (references 
in Barnes 2003). In the mass range concerned by the AO-FO 
stellar sample of Abt & Morrell (1995), the detected excess of 
low rotators was assigned entirely to Am and Ap stars. The Ap 
phenomenon is, however, evolutionary related, because it ap- 
pears after the stars completed at least the first third of their 
MS life (Hubrig et al. 2000; St?pieh 2000). On the other hand, 
most Am stars do not have luminosity class determination, so 
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Fig. 10. Distributions of angular velocities £l/Slc derived from v distributions in Fig. 7 and using Appendix D. Only the variability band of each 
PDF is displayed. 



that their evolutionary status is unclear. Moreover, a great frac- 
tion of them, if not all of them, are binjiries (Debemardi 2000). 
The Am+Ap group cannot then be considered witnessing a 
neat low velocity component in a possible bimodal distribu- 
tion. We cleansed our stellar sample from objects with known 
Am and Ap signatures, and from all known close binaries to 
avoid strong spurious effects on the velocity distributions car- 
ried by braking effects that can not be considered forming part 
of the initial MS conditions for single stars. Thus, as binaries 
and stars with low rotation possibly due to evolutionary effects 
were removed, the distributions we found for field B9, AO-Al 
and F0-F2 stars may be considered as genuinely bimodal. 

The simplest explanation that can be tempted to these bi- 
modality is that the braking of surface rotational velocities 
trough interactions with the circumstellar disc prior to the 
ZAMS phase is more or less lasting. The two velocity peaks 
would then imply the existence of two preferential interaction 
time scales, which is difficult to prove. However, this explana- 



tion sounds not sufficient, because in stars with spectral types 
from A2 to A9 the bimodal aspect of velocity distribution is 
not obvious. Perhaps, also some phenomena of internal angu- 
lar momentum redistribution can be present. 

It is worth noting the lack of stars with n/Qc ^ 0.4 in the 
A4-A6 spectral-type group. This group is roughly in the mid- 
dle of the transition from stars with radiative envelopes to those 
with convective envelopes. The behavior of the 1600 A inten- 
sity jump (Bohm-Vitense 1982) in stars from Al to A7, reveals 
a strong ionization balance change in the envelope, and so in 
their structural patterns, as we pass from one spectral type to 
another. It can thus happen that either the PMS angular mo- 
mentum setting-up phenomena are somehow related with stel- 
lar structure and its external ionization balance, or that in the 
early ZAMS phases some internal redistribution of angular mo- 
mentum operates making that f2/Qc ^0.4 is unliliely. 

On the other hand, in all subgroups there is a rough aver- 
age of 15 % stars with Q/Oc ^ 0-8. These stars are subjected 
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to rotationally induced geometrical distortions and to gravita- 
tional darkening. Since the local effective temperature is, ac- 
cording to the von Zeipel theorem, T^g oc g (g is the local 
stellar surface gravity) in stars with radiative envelopes and 
T^gccg^-^ (Lucy 1967; Perez Hernandez et al. 1999) in objects 
with convective envelopes, a strong mixing of stars with differ- 
ent masses is expected in the hottest stellar groups. The grav- 
itational darkening effect can also play a role in the selection 
of objects in the A4-A6 spectral-type group. In fact, pole-on 
stars rotating at rates Q/Qc ~ 0.8 which are seen as A3 (aver- 
age hemisphere effective temperature T^^'^ ^ 8700 K), would 
be seen equator-on likely as A7 (average hemisphere effective 
temperature T^^"^'" ^ 7900 K). In non-rotating stellar models, 
the effective temperature range 7900 ^ T^g 8700 K corre- 
sponds to the transition from stars with radiative envelopes to 
those with convective envelopes. It would then be interesting to 
investigate if such objects present a radiative/convective struc- 
tural dichotomy and if the lack of low rotators in the A4-A6 
group is somehow related with this structural characteristic. 

Since there is a correlation between the strength of mag- 
netic fields, low rotation and presence of the CP character in 
A to B-type stars (Mathys 2004), a thorough inquire on the 
nature of objects with low rotation in our sample that do not 
have apparent CP aspects appears timely to address. However, 
the separation of stars in our sample with genuine low equato- 
rial velocity v from those with small sin / is not obvious. Our 
planned determination of detailed fundamental parameters and 
identification of fast rotating stars with low sin i with model at- 
mospheres for fast rotators may perhaps be of some help to this 
purpose. 

6. Summary and conclusion 

This work was based on a large and homogeneous sample of 
vsin/ measurements of B9 to F2-type gathered from Royer 
et al. (2002b) and Abt et al. (2002). This sample was carefully 
cleaned of evolved stars and known chemically peculiar and/or 
binary stars, resulting in a sample of about 1 100 main-sequence 
single stars. This stellar sample was divided into six spectral 
type groups, chosen so as to warrant in each statistical signif- 
icance and mass resolution, to allow for detection of possible 
mass dependencies of velocity distributions. 

The distributions of true equatorial velocities have been 
obtained from observed v sin / distributions using the iterative 
technique described by Lucy (1974) to rectify from the error 
and random inclination distributions. 

- The mean true equatorial linear rotational velocity of A2 to 
F2 MS dwarf stars are higher than quoted in the literature. 
The differences with the previous estimates can attain Av ^ 
50kms 

- The distribution of rotational velocities of B9 to Al-type 
"normal" field dwarf stars are genuinely bimodal. This find- 
ing is contrary to the results obtained by Guthrie (1982), 
Abt & Morrell (1995) and Abt (2000). The presence of 
modes around 50 and 200km s"' may be due to formation 
processes and phenomena of angular momentum loss and 
redistribution undergone during PMS phases. 



- There is a striking lack of rotators with v ^ 50-70 km s"' 
from A2 to A9-type MS field stars. The minimum of stars 
in this interval is for A4-A6-type stars. 

- The F0-F2 group seems to have a bimodal distribution. 

- A great fraction of A-type stars are fast rotators: Q/Qc ^ 
0.5. Hence, it should not be excluded that the faster rotating 
A-type stars, mainly A4-A6-types, might have envelopes 
that have non homogeneous structures in latitude: radiative 
(polar regions) and convective (equatorial regions). 

Further detailed analysis of the present stellar sample using 
fundamental parameter determinations and models of rotating 
stars will be given in forthcoming papers. 
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Appendix A: Notes on common v sin i standard 
stars 

The four conmion B-type stars between measurements from 
Papers I and 11 and the standard values from Slettebak et al. 
(1975) are listed in Table 1. There are few available recent data 
about any v sin i measurement for these objects. Some details 
are given for them below: 

- yS Tau has been observed at high angular resolution by 
Richichi & Percheron (2002) and no suspicion of binarity 
has been found. The resulting uniform disk angular diame- 
ter is 1 .56 ± 0. 1 1 mas. 

- T Her was discovered by HIPPARCOS to be a slowly pul- 
sating B star with a period of l'?24970 ± 0?00008 (ESA 
1997). Masuda & Hirata (2000) recently observed its line- 
profile variation on He 1 447 1 and Mg n 448 1 A. These vari- 
ations could induce some measurement effects on the v sin ; 
and explain the discrepancy found between this work and 
the quoted results in the literature. 

- y Lyr is described in Appendix of Paper II. 

- (j/ Aqr is described in Appendix of Paper I. 

Appendix B: Error on the empirical dispersion 

The empirical variance is estimated by = 2; where 
X, are supposed to have a dispersion cr. The variance of this 
estimator is: 



var(i^) = n ^ ^ var(x^). 



The dispersion of x? is 

(T^2 ss 2 \Xi\ cr, 
then var(x^) x4xfcr^. 



(B.l) 



(B.2) 
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It easily follows that 



(B.3) 



because lim„^oo s = cr. Thus 0-^2 -2cP- iT^I^. 

The dispersion of the empirical standard deviation s is 
therefore: 



a, « cr,2(2i) 



an oc cr. 



(B.4) 



Q.e.d. 



Appendix C: Notes on manipulating distributions 

C.1. Kernel method 

This method for density estimation is extensively described 
in Silverman (1986) and Bowman & Azzalini (1997). 
Applications can also be found in the literature (see e.g. Arenou 
1993; Jorissen et al. 2001). Using the kernel function K, the 
probabihty density function O is estimated by: 



1 -A 1 jx-Xi\ 

1=1 



(C.1) 



where n is the sample size and h is the band width parameter. 
The mean integrated squared error (MISE) is used to measure 
how efficient O is in estimating O. It is defined as: 



MISE((t>) = £ (D(x) - <l)(x) dx, 



(C.2) 



where E{f{x)) denotes the expectation value of a function /(x) 
in a variable x. Minimizing the MISE allows the choice of the 
best kernel function. Although it is well known that the so- 
called Epanechnikov kerneP is the kernel function minimizing 
the MISE (Silverman 1986), the Gaussian kernel, with a lower 
efficiency, has been chosen in order to make use of the Sheather 
& Jones (1991) method for estimating the bandwidth parameter 
h. This latter parameter is estimated using the Fortran subrou- 
tine available from Sheather & Jones (1991). The kernel es- 
timator O is computed using the C code written by Arenou 
(1993). 

The variability bands, as described by Bowman & Azzalini 
(1997, Chap 2.3), are computed to assess the significance 
of the modes in the derived distributions. For a given esti- 
mated density function O defined by Eq. C.l, and according to 
Bowman & Azzahni (1997), the variabihty band width around 
V5is: 



2e = 



nh)-^ J K\i)At. 



(C.3) 



^ The Epanechnikov kernel is defined by: 
\ otherwise. 



In the case of a Gaussian kernel function (Eq. 5), it becomes: 

2e = ^^(nftr^~(2A^ 
. . 0-2656 

therefore e « — (C.4) 



This is applied to the estimated PDFs <l)iog(ln(v sin /)) in 
Sect. 4.1, and the corresponding e values are given in Table 5 

for each estimated ^^iog. The variability band of Oiog is 

bounded by both functions ^ -^j^iog ± f j ■ 

These bands, first computed for the smoothed estimated 
distributions, are propagated through the whole processing 
(Lucy-iteration technique) from smoothed density functions to 
final deconvolved ones. 



C.2. Lucy iterative teciinique 

Lucy (1974) gives an iterative deconvolution scheme to solve 
the problems that consist in estimating the density function 
giY) of a quantity Y from a sampled population characterized 
by: 



/(x)= g(y)Px{x\y)dy, 



(C.5) 



The first guess of g(y), g'^iy), is taken to be the uniform law. 
According to Lucy, the iterated r-th functions are: 



(y)Pxix\y)dy, 



rn 
J f' 



fix) 
fix) 



Px(x\y)dx, 



(C.6) 



(C.7) 



where f(x) is the estimation of the PDF f(x) based on the ob- 
served data. 

These operations are carried on using the smoothed distri- 
butions and overlapping parabolas integration method (Davis 
& Rabinowitz 1984). 

To stop the iterations, Lucy used a tsst. We pre- 
ferred to use a Kohnogorov-Smimov test to check whether 
the hypothesis /'' = f is true within a 1 %-threshold. In 
the calculations carried out on the vsinj data and detailed 
in Sect. 4.1, the number of iterations of the Lucy technique 
is 3 or 4 for the deconvolution by the errors distribution 
[/(x) = <l)iog(hi(vsinO),g(>') = ^iog(hi(i?))], and ranges from 
12 to 40 for deconvolution by the inclinations distribution 
[fix) = ^i-»),g(y) = riv)l 

Appendix D: Angular velocities Q/Qc 

In this representation, the variation of the polar radius with cj = 
Q/Qc is neglected. 

The Roche model of surface equipotential for a rotating star 

gives: 



Rc 



1 + 0.5 s, 



(D.l) 
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where is the critical radius: Rc = Re(<o = 1), and 




(D.2) 



Asw - o) it follows that: 

a) = w[l-0.5w^). (D.3) 

proxima.obspm.fr:/home/royer/ScienceA^sini/in/astro-pli^, 
more vsinilllcorB.bbl 
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